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Abstract. Building on the work of lSchweizerl (|1995l ) and lCernv and Kallsenl (|2007l ) 
we present discrete time formulas minimizing the mean square hedging error for 
multidimensional assets. In particular, we give explicit formulas when a regime- 
switching random walk or a GARCH-type process is utilized to model the returns. 
Monte Carlo simulations are used to compare the optimal and delta hedging meth- 
ods. 

Hedging, option pricing, GARCH, regime-switching 



1. Introduction 

In many applications, one is interested in finding a discrete-time dynamically traded 
portfolio so that its value at maturity is as close as possible to a target function of 
the underlying assets. When the target function is a payoff, this can be interpreted 
as option pricing and hedging. However, sometimes the target function is not a 
payoff as it happens when one tries to replicate hedge funds or create synthetic funds 



( IPapageorgiou et al 



2008). As suming that the error measure is the average quadratic 



hedging error, ISchweizerl ( 119951 ) solved the hedging problem for one risky asset. He 
showed that the initial value of the portfolio, which can be interpreted as the "value" 
of the option, is the average, under the "real probability measure", of the discounted 
payoff, multiplied by a martingale, which is not necessarily positive. In the latter case, 
the martingale cannot be used as the density of an equivalent martingale measure. 



Partial funding in support of this work was provided by the Natural Sciences and Engineering Re- 
search Council of Canada, by the Fonds Quebecois de Recherche sur la Nature et les Technologies, the 
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Even if the hedging problem has been solved quite generally by 



the o ne- dimensional case , it seem s to have been ignored or forgotten, e.g., (IBouchaud and Potters! . 



2002|) or 



Schweizerl (|19951) in 



Cornalba et al. 



( 120021 ). More troubling, delta hedging, based on the Black- 
Scholes-Merton model, is still used in practice even if it has been shown that the 
geometric Brown i an m otion model is an inadequate model for the underlying assets 



( iKat and Palaro 



20051 ) . Even when the geometric Brownian motion model is ade- 



quate, the hedging error in discrete time is not zero, co nverging only to zero as th e 



number 


0: 


' hedg 


Wilmott 




2006, 



hedg ing periods tends to infinity. See, e.g. 



Boyle and Emanuel! ( I1980l ). 



Chapters 46-47). In addition, when the market is not complete 
but there is no arbitrage, there are infinitely many equivalent martingale measures 
(EMM). One then has to choose the "best" martingale measure with respect to some 
utility criterion. There is a huge literature o n this subject . Indeed, when a NGARCH 



process is utilized to model the log-returns, 



Duanl (119951 ) proposes a solution to the 



EMM problem. Unfortunately, Duan also suggests a delt a hedging strategy, which 



has been shown to be wrong by 



Garcia and Renault 



(11998 ). 



Papageorgiou et al. 



(120081 ) 



Motivated by applications in hedge fund replication 
proposed a locally optimal solution minimizing the average quadratic hedging er- 
ror at each period for the general multidimensional asset case. They erroneously 
claimed that it was globally optimal, which is only true if the discounted underlying 
assets are martingales. A first motivation for the p resent paper is to correct that 



mistake and give ex plicit formu l as fo r the results in 
generalizing those of 



Cerny and Kallsenl ( 120071 ) and 



Schweizerl ( 119951 ). A second motivation is to show that, when 
regime-switching random walks and GARCH processes are used to model the returns, 
the optimal solution of the hedging problem yields superior results to those obtained 
by delta hedging. 

The optimal solution of the discrete time h edging problem is d e scribe d in Section 



121 giving explicit expressions for the results of 



Cerny and Kallsenl ( 120071 ). It is worth 
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noting that when the asset value process is Markovian, or a component of a Markov 
process, the optimal solution can be implemented using approximation techniques 
of dynamic programming. Two such cases are considered. Finally, in Section [31 
simulations are used to compare optimal hedging with delta hedging for geometric 
random walks and NGARCH models. 



2. Optimal hedging strategy in discrete time 

Denote the price process by S, i.e., Sk is the value of the d underlying assets at 
period k and let F = {J~k, k = 0, . . . , n} be a filtration under which S is adapted. 
Assume that S is square integrable. Set A k = fikSk — Pk-iSk-i, where the discounting 
factors /3k are predictable, i.e. /3k is J-fc_i-measurable for k = 1, . . . ,n. 

The aim of this section is to find an initial investment amount Vo and a pre- 
dictable investment strategy (f> = (<f>k)k=i sucn that T is square integrable and 



which minimizes the expected quadratic hedging error E 



G = G Vr 



G(y , 



, where 



(3 n C - V n , and V k = V + Eli 4>j^j, k = 0, 



n. 



The exi s tence of an optima 



Schweizerl (119951 ). 



Motoczyhski 



soluti on was proven first in the univariate case by 



nishing explicit solutions. Finally, 



(120001) considered a mu l tivaria te setting without fur 



Cerny and Kailsenl (120071 ) treated a much more 



general case. However, in the discrete time case, it is faster and easier to find directly 
the explicit formulas than trying to recover them from their results. 

2.1. Offline computations. Once a dynamic model is chosen for the asset prices, 
one must start with some computations that are necessary for the implementation. 
Set P n+1 = 1, 7„, + i = 1, and for k = n,...,l, define A k = E (A k AjP k+1 \ T k -i) , 
lik = E (AfePfe+ilTfc-i), b k = A/: 1 ^, P k = rijLfc (l and 7 fc = E(P k \J r k^i), 

provided these expressions exist. Under some extra assumptions given below, it can 
be shown that they are indeed well defined. The proof is given in Appendix IA.U 
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Lemma 2.1.1. Suppose that E(pj k+ i\F k _i)A k — fikfij is invertible P-a.s., for every 
k — n, . . . , 1. Then 7& G (0, 1] and A k is invertible for all k = 1, . . . , n. In addition 
(7fc + i)£ =0 is a positive submartingale. 



Remark 2.1.2. In the univariate case, 



Schweizen u 99a) states sufficient conditions 



for the validity of the assumptions of Lemma \2.1.1[ It is not obvious how they could 
be generalized to the multivariate case. Therefore, in most applications, one has to 
verify these conditions, often using brute force calculations. 

2.2. Optimal solution of the hedging problem. 

Theorem 2.2.1. Under the assumptions of Lemma \2.1.1\ the solution ^Vo, 0^ of 
the minimization problem is Vq = E{fi n CPi)/ r )i, and <p k = a k — Vk-i°k, where a k = 
A~ k x E (f3 n CA k P k+1 \ F k -i), k = n,...,l. 

The proof is given in Appendix IA.21 

2.2.1. Option value. Let C k be the optimal investment at period k, so that the value 
of the portfolio at period n is as close as possible to C, in terms of mean square 
error. One could then interpret C k as the value of the option at period k, for any 
k = 0, . . . , n. It then follows from Theorem 12.2.11 that C k is given by 

(2-2.1) fi k C k = — — — — , k = 0,...,n, 

-CH-Tfc+lK k) 

so one can write 

(2.2.2) = —E {(3 k C k (1 - bj A fc ) 7fc+i|^-i) = E {(3 n CU k ■ ■ -U n \7 k . x ) , 

7fc 

where U k = jJpjr^~~j > k = 1, . . . ,n + 1, while an alternative expression for a k is 

(2.2.3) a k = A~ l E {f3 k C k A klk+l \T k ^) . 
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Remark 2.2.2. Setting Z = 1 and Z k = YLj=ity> ^ = 1 , - - • , ^ ? one obtains that 
(Z k , (3 k C k Z k , (3 k S k Z k )^ =0 are martingales. However, in most applications, Z does not 
define a change of measure unless it takes only positive values. 



2.2.2. Implementation issues. If the process S is Markov and C n = C n (S n ), then 
Ck = Ck(Sk), ctk = ctk(Sk-i), and b k = b k (S k -i). It foll ows that all these function s 



Papageorgiou et al 



fl2008h 



can be approximated using the methodology developed in 
Another interesting case is when S k is not a Markov process but (Sk, h k ) is, even if 
h k is not observable, as in GARCH and regime-switching models. In this case, if 
C n = C n (S n ), then C k = C k (S k ,h k ), a k = a k (S k -i,h k _x), b k = b k (S k _x, and . 

Ik = 7fc(<Sfc-i, hk-i), for k = n + 1. All these functi ons can be approximated 



using the methodology developed in 



Remillard et al. 



(120101 ) for the regime- switching 



case. Implementation of the hedging strategy then requires predicting h t . 

Remark 2.2.3. One could suggest to use the smallest filtration to get rid of the 
unobservable process h but in this case, all conditional expectations based on T k would 
depend on all past values S , . . . , S k , making it impossible to implement in practice. 



2.3. Verification of the assumptions of Lemma 12.1.11 In what follows, we con- 
sider some interesting models used in practice, for which it is possible to show that the 
assumptions of Lemma 12 . 1 . II hold true and that the optimal solution can be computed 
via a dynamic program. 



2.3.1. Regime-switching geometric random walks. An interesting model, which in- 
cludes geometric random walk models, is to consider a regime-switching geometric 
random walk. Theses models can display serial dependence in the log-returns and 
may account for changing volatility over time . For implementation issues, including 



estimation and goodness-of-fit tests, see, e.g., 



Remillard et al. 



(120101 ). 
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To define the process, suppose that r is a finite homogeneous Markov chain with 
transition matrix Q with values in {1, . . . , / } representing the non-observable regimes 
and set (3 k S l k = S l JlLii 1 + $}) i = 1, • • • , where, given n = z'i, . . . , r n = i n , 
£i, . . . , £ n are independent with ^ ~ Pj., j = 1, . . . , n, Ei(£j) = (Cj I T i = i) = m(i), 
and Ej (Cj&f") = B{i). The interpretation of the model is easy: At a given period 
t, a regime r t is chosen at random, according to the Markov chain model, and given 
Tt = i, £t is chosen at random according to distribution Pj. When there is only one 
regime, one obtains a geometric random walk where all £s are independent. 

We assume that the B(i) — m(i)m(i) T is invertible for any i = Setting 
X k = (3 k S k , one gets A k = X k - X k -i = D(X k _i)£ k , k = l,...,n, where D(s) is the 
diagonal matrix constructed from vector s. Note that S is not a Markov process in 
general but (S, r) is a Markov process. The validity of the assumptions of Lemma 
12. 1.11 follows from the next result, proved in Appendix I A. 31 



Proposition 2.3.1. For any k = 1, . . . ,n and i = 1, . . .,1, ^ k = 7^(^-1), 7fc(i) £ 
(0, 1], A k = A k (S k -i,T k -i) and b k = b k (S k -i, r k -i), where 

1 



(2.3.1) 



A k (s,i) = fi_ x D{s) \J2Qijlk+i(j)B(j) \ D(s), 

.i=i 



(2.3.2) 



b k (s,i) = D l ((3 k _ 1 s)p k+1 (i), 



(2.3.3) 



7fcW 



^Qijlk+iU) {l - Pk+i(i) T m{])} 



with p k +i(i) = {EJ=i Qijlk+iti)B(])} {EJ=i <5ij7fe+i(j) m (j)}- 
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If in addition C = C(S n ), then C k = C k (S k ,r k ) and a k = a k (S k -i, Tfc-i), where 



Cfc-l(s,z) = -q-^-J^Q 



x / C fc 



^D(s)(l + y), A {1 - Pfe+1 (z) T y} P,(dy), 



I I" 1 I 



k-l 



. i=i J j=i 

2.3.2. GARCH-type models. Here, one assumes that A fc = (3 k S k —f3 k -S k -i = /3 k -iS k -i£ k , 
with £k = 7r 1 (/ifc_ 1 , ejt), and ft^ = ^(/i^-i, e&) with 7r2 having values in some set and 
where the innovations e k are independent and identically distributed with probability 
law v. It is immediate that (S k , h k ) is a Markov process. Furthermore, almost all 
known GARCH(1,1) models can be written in that way. 

Suppose that for every given possible h E?i, ni(h,y) is not constant i^-a.s. Using 
Proposition IB. 0.11 and reverse induction, as in the proof of Proposition 12.3.11 it is 
easy to show that the assumptions of Lemma 12.1.11 are met, and that for all k = 
n, ...,1, 7fc = 7fc(^fc-i) and A k (s,h) = /3l_ v s 2 B k (h), b k (s,h) = j^^j^ , where 
B k(h) = $ -K{(h,y)~j k+l {K 2 {h,y)}v(dy), m k {h) = J m(h, y)lk+i {ir 2 (h,y)} v(dy), and 
lk(h) = J |l - 1[ §^ni{h,y)\-i k+1 {Ti2{h,y)}u{dy). Also, if C = C n (S n ), then 



C k „x(s,h) = I C k 



§ 

—{l + 7n(h,y)},7r2(h,y) 
Pi 



X7fc+i {n 2 (h,y)} { 1 - ^^-^(h, y) \ v{dy), 
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g 

— {1 +n 1 (h,y)},ir 2 (h,y) 
Pi 



xik+i {n 2 (h, y)} 7Ti(/i, y)v{dy). 
Hence, the optimal solution can be written as a dynamic program. 

3. Examples of application 

In this section we consider pricing and hedging of European calls for two geometric 
random walk models, when the returns are i.i.d. Gaussian and i.i.d. differences of 
Laplace distributions, and for a NGARCH model. It follows from the previous sections 
that optimal hedging solutions exist for these cases, and the optimal solution can be 
written as a dynamic program associated with functions of a finite number of variables. 
For solving such dynamic programs, we discretize the state space into a finite grid and 
we compute approximations of expectations using Monte Carlo simulations at every 
point of the grid. Linear interpolations are used for points outside the grid at each 
time step. Since the expectations are always with respect to the same probability 
measure, only one sequence of random numbers may be used, using the ideas in 



Del Moral et al. 



(12006 



2012)) 



3.1. Geometric random walk models. Here we consider discretized versions of 
the Black-Scholes (BS) and Variance Gamma (VG) models for the underlying asset 
over 23 periods. In each case, the 22 periodic returns are i.i.d., so that the mean and 
volatility at maturity are respectively 9% and 6%. For the BS model, the returns 
are Gaussian, while for the VG model, the returns are differences of i.i.d. Gamma 
variates, so that the distribution at maturity is Laplace (double exponential). These 
models are particular cases of regime-switching models with only one regime. We 



do no t consider regime- switching models since it has been done in 



Remillard et al. 



( 120101 ) . where the daily log-returns of the S&P 500 are analyzed. 
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We are going to price and hedge a call with strike K = 100 and maturity 1 year, 
using 22 replication periods and a 2000 points discretization of the asset values over 
the interval [80,120]. The annual rate is 5%. A sequence of 10000 random points 
were used for the computation of functions and Ajt, while 10000 paths were used 
to compute the hedging errors. Delta hedging is optimal in the continuous time limit 
for the BS model, but not for the VG model. As expected, according to Figures (TH2J 
the values of the call Co and initial investment strategy <fti, obtained from the optimal 
hedging, are close to those obtained using the Black-Scholes formula (even with 22 
hedging periods), while they differ for the VG model. This is also reflected in the 
distribution of the hedging errors, as illustrated in Figure O 




Q l -~ i 1 1 1 1 1 1 -5 1 1 1 1 1 1 1 1 1 

80 85 90 95 100 105 110 115 1 20 80 85 90 95 100 105 110 115 1 20 



Figure 1. Call option value Cq for the Black-Scholes (left panel) and 
Variance Gamma models (right panel) with 22 periods of hedging. 



Initial investment strategy tor Black-Scholes model Initial investment strategy for Variance Gamma model 




80 85 90 95 1 00 105 110 115 1 20 80 85 90 95 100 105 110 115 1 20 



Figure 2. Initial investment strategy 0i in the underlying asset for the 
Black-Scholes (left panel) and Variance Gamma models (right panel) 
with 22 periods of hedging. 
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Density estimation of the optimal hedging errors: Payoff-Portfolio 




■ Optimal hedgind 
- Delta-hedging ] 



Density estimation of the optimal hedging errors: Payoff-Portfolio 



Optimal hedgind 

Delta- hedging ] 



Figure 3. Estimated densities of the hedging error G for the Black- 
Scholes (left panel) and Variance Gamma models (right panel) with 22 
periods of hedging. 

Descriptive statistics of the hedging errors are given in Table [TJ Simulations can 

also be used to show that as the number of hedging periods increases, the hedging 

error tends to zero for the BS model, while it is never for the VG model. Note that 

the RMSE of the optimal hedging is always less than the one of the delta hedging. 

Table 1 . Statistics of hedging errors (Payoff-Portfolio) for the Black- 
Scholes and Variance Gamma models. 





Black-Scoles 


Variance 


Gamma 


Stats 


Optimal 


Delta 


Optimal 


Delta 


Average 


-0.0065 


0.0076 


0.0151 


0.0518 


Median 


-0.0014 


0.0029 


-0.1543 


-0.2886 


Volatility 


0.2537 


0.2774 


1.0510 


1.2529 


Skewness 


1.3781 


0.5515 


6.1479 


6.0965 


Kurtosis 


22.0975 


8.8391 


63.3897 


60.9747 


Minimum 


-1.7986 


-1.7578 


-4.4372 


-1.5160 


Maximum 


3.5389 


2.2978 


18.5163 


23.6422 


VaR(99%) 


0.8231 


0.9244 


4.6976 


5.9450 


VaR(99.9%) 


1.9576 


1.5779 


11.7905 


12.2222 


RMSE 


0.2538 


0.2775 


1.0511 


1.2540 



3.2. NGARCH model. As in 



Duanl (119951 ). we consider the NGARCH model 



where e €fe - 1 = r + \^Jh k -i - \h k -i + y/h k -ie k , and h k = a + atih k - X e 2 k + ft/ifc-i, 
with e k ~ N(0,1) and parameters a = 1.524 x 10" 5 , a x = 0.1883, ft = 0.7162 
and A = 7.452 x 10" 3 . Under the EMM, we have e ik - 1 = r + + ^/h k _ x ^k 

and h k — ao + a\h k _i(e k — A) 2 + j3ih k -i- We price and hedge a call with strike 
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K = 100 and maturity 30 days using daily replication, using a grid of 500 points 
for the asset on [60, 140], while the grid for the volatility consists in 90 points of the 
interval [.00005, .0007]. The annual rate is 0%. In what follows, B&S hedging means 
delta hedging using the B&S formulas, while Duan's methods consists in picking his 
suggested EMM and taking the delta of the option. The value of the option and 
the initial number of asset are displayed in Figure HJ while descriptive statistics of 
the 10000 hedging errors are given in Table [2] for the three hedging methodologies, 
showing that the errors are more concentrated about for the optimal hedging. 



Call option value for NGARCH model lnitial investment strategy ( $ ^ 




Figure 4. Optimal hedging call option value Co and initial investment 
strategy <pi for the NGARCH model with 30 periods. 



Table 2. Statistics of hedging errors (Payoff-Portfolio) for the 
NGARCH model. 



Stats 


Optimal 


Delta 


Duan 


Average 


-0.0159 


-0.0954 


0.0085 


Median 


-0.1549 


-0.2094 


-0.1297 


Volatility 


0.8568 


0.8951 


0.9032 


Skewness 


1.7205 


2.7558 


2.9171 


Kurtosis 


10.5790 


29.3709 


29.9947 


Minimum 


-1.9966 


-2.2302 


-2.0789 


Maximum 


9.9114 


17.4873 


17.7232 


VaR(99%) 


2.7698 


2.8613 


3.0831 


VaR(99.9%) 


5.4893 


6.3574 


6.5070 


RMSE 


0.8569 


0.9001 


0.9033 
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4. Conclusion 

In this paper we presented the optimal solution for a discrete time hedging port- 
folio. When the underlying process is Markov or a component of a Markov process, 
the optimal hedging strategy depends on deterministic functions that can be approx- 
imated. We also find explicit formulas for two interesting models. Finally, numerical 
simulations show that optimal hedging is preferable to delta hedging. 
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Appendix A. Proofs of the main results 

A.l. Proof of Lemma l2.1.1l First, we will show that 7^ G (0, 1] and Ay. is invertible 
for all k — 1, . . . , n. By hypothesis, E(jk+i\ J^k-i)^ — Hk^l is invertible for all k = 
1, . . . , n. In particular, it is true for k = n, yielding that S n = A n — fi n fi^ is invertible, 
which is the conditional covariance matrix of A n given T n -\. It then follows from 
Proposition IB. 0.11 that A n is invertible. Without loss of generality, one may assume 
that A n is diagonal. Otherwise, we diagonalize it in the form A n = M n B n Mj , with 
M n ,B n J^.j-measurable, B n is diagonal, M n = I and set A n = MjA n . Since 
M n Mj = I, it follows that M n is bounded, so A n is square integrable. Finally B n = 
E(A n A^\ J s n -i). A n being diagonal, it then follows that b^A n is square integrable 
and 7 n = 1 — b^fj, n = 1 — fj^A^ 1 ^. It also follows from Proposition IB. 0.11 that 
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.Sl^l » SO 7n = 1+u tU u £ (0, !]■ AS a reSult 5 7n < ^(7n+l^n) = 1. 

The rest of the proof follows easily by reverse induction, using Proposition IB.O. II with 
the mean and covariance matrix of A k under the probability distribution Q k , with 
Q k (0) = E(I j k +i\J 7 k -i)/E(j k+l \J : k ^ 1 ), O e T u for k = n - 1, . . . , 1. □ 



A. 2. Proof of Theorem 12.2.11 Using the proof of Lemma I2.1.1[ one can easily 
check that a k , b k and <p k make sense and that <flA k is square integrable. Next, it 
is easy to check that a necessary and sufficient condition for (Vo, j to minimize 



E 



G Vr 



all k 



is that E |G (Vo, ~$\ } = and E \g (V , ~$\ A fc |.F fc _i} = for 
n. The necessity comes from the fact that for any event O G J-" k -i, 



one must have 



— I W 

de \e=0 



'G(y , <f>) -eI G A fc } = -2e{g(v , </>) A fc I G }, 
which is equivalent to the condition E |(jr (Vo, 4> \ A k \ j = 0, while the condition 
E |g (Vo, 0j|=O comes from the fact that for any 6, one must have 



de 



E 



e=0 



2E<G[V( 



To see that the conditions are sufficient, it suffices to check that 



E 



[G(y, + e^JTpj] 2 =e {G(y ,~$)Y + £|U + X>Ja 




The proof that is the solution is based on the following equation, which can be 
easily proven by induction. 



(A.2.1) E(V n \F k ) = V k E{P k+l \F k ) + E {/3 n C(l - P k+1 )\F k ) , k = l,...,n. 



To complete the proof of theorem, note that from (1A.2.1I) , 



(A.2.2) E[G{y Q ^) \T k } = E{p n CP k+1 \JF k )-V k E{P k+1 \JF k ), k = 0,...,n. 
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Using fM2| , one has E {g((v , ^) A fc | T k } = E^ n CA k P k+l \T k )-E(V k A k P k+1 \T k ), 
so E{G((y ,^ Ak\T k -i) = E(f3 n CA k P k+1 \T k ^) - E(V k A k P k+1 \F k ^) = A k {a k - 
V k -!b k - <p k ) = 0. Hence E {g((v , ~P) } = E(P n CP 1 ) - V E(P 1 ) = 0. □ 

A. 3. Proof of Proposition 12.3.11 The result is obviously true for k = n + 1. 
Suppose it is true for k + 1. For i given, set iij = <3ij7fc+i(j) / D, where D = 
X)j=i Qijlk+i(j)- By hypothesis, 7Ti, . . . , 717 are probabilities adding to 1, so if X ~ Pj 
with probability 7Tj, then •y k (i) = D (l — n T B~ l ji} , where /i = -E(AT) and -B = 
E (XX T ) . Let £ be the covariance matrix of X. It is non singular since the co- 
variance of X under Fj is non singular. It then follows from Proposition IB. 0.11 that 
1 — /i T £> _1 /i = 1+ t^-i^ > 0. Since D > by hypothesis, one may conclude that 
7fc(i) > 0. As a by-product we get that ^y k (i) < 1 if 7fc+i(j) < 1 for all j = 1, . . .. 
Since that is true for 7 n+1 = 1, one may conclude that for all k — 1, . . . , n, 7fe(i) < 1. 
The rest of the proof is easy. □ 

Appendix B. Auxiliary results 

Proposition B.0.1. Suppose A = S+6fe T where £ is symmetric and invertible. Then 
A is invertible, and A' 1 = - ^rpr^ ■ Moreover, 1 - b T A^b = ^ i 1 ^^ > 0. 

Proof: Since A (V 1 - X+b%-^b ) =I,Ais invertible and A" 1 = E" 1 - ^re- 
setting c = b T ^-% one gets 1 - b T A~ l b = l- c+ I ^ = T ^>0. □ 
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Initial investment strategy for Black-Scholes model 
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Initial investment strategy for Variance Gamma model 
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